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We show that the time development of the Coulomb gap in a Coulomb glass can involve very 
long relaxation times due to electron rearrangement and hopping. We find that an applied magnetic 
field reduces the rate of electron hopping and, hence, Coulomb gap formation. These results are 
consistent with recent conductance experiments on thin semiconducting and metallic films. 



The competition between interactions and disorder re- 
sult in glassy dynamics that are often associated with 
very long relaxation times extending over many decades. 
One might not expect the same to be true in an electronic 
system since electrons typically respond very quickly. 
However, in this paper we show that in the presence of 
strong disorder, electrons can indeed have very long re- 
laxation times. This occurs in a Coulomb glass which is 
an insulator with randomly placed electrons that have 
Coulomb interactions. Heavily doped semiconductors 
and disordered metals are examples of such systems. 
Coulomb interactions between localized electrons result 
in a so-called Coulomb gap in the single particle density 
of states that is centered at the Fermi energy [|]-|| . We 
have done a calculation in which we follow the time de- 
velopment of the Coulomb gap. In order to produce this 
gap, electron rearrangement must occur and the associ- 
ated hopping can involve very long time scales. 

These long relaxation times are consistent with recent 
experiments on thin semiconducting [^|J|] and metallic 
1^ films which have shown that in the presence of strong 
disorder, electronic systems can relax very slowly. These 
films were grown on insulating substrates which sepa- 
rated them from a gate electrode that regulated the elec- 
tron density, and hence the chemical potential, of the 
film. The conductance G was measured as a function 
of the gate voltage Vc- If Vq sat at a particular value, 
Vo, for a long time and then was varied over a range of 
voltages, there was a dip in the conductance centered at 
Vo . We identify this dip with the Coulomb gap in the 
density of states because the value of the conductance de- 
pends on the density of states at the Fermi energy . 
In Mott's picture of variable range hopping, the hopping 
conductivity increases when the density of states at the 
Fermi energy increases, since there are then more states 
to which an electron at the Fermi energy can hop 
We identify sweeping Vq with varying the chemical po- 
tential without allowing time for equilibration. In effect 
the sweeps scan the density of states. Thus we expect 
the conductance to increase with the density of states 
and hence as the gate voltage Vg moves away from Vo- 

In the experiments, if the gate voltage was changed 
suddenly from, say, Vo to Vi , the conductance had a very 



fast initial rise, followed by a period of rapid relaxation, 
which in turn was followed by a long period of very slow 
relaxation. In some cases the relaxation was logarithmic 
in time. Our interpretation of this is that when the gate 
voltage is changed, the Fermi energy changes, and time 
dependent relaxations arise because the system must dig 
a new hole in the density of states at the new Fermi 
energy and remove the old hole at the old Fermi energy. 
Indeed, Vaknin, Ovadyahu and PoUak j|] found that sub- 
sequent sweeps of the gate voltage revealed that the old 
dip in the conductance at Vo fades with time while a new 
dip centered at Vi increased with time. The dip in the 
conductance and the long time relaxation were present 
only at very low temperatures, not at higher tempera- 
tures (T ~ 20 K). 

Ovadyahu and PoUak also found that in a magnetic 
field {H = 9 Tesla) the long time relaxation rate associ- 
ated with a change in the gate voltage was even slower. In 
addition to the spin mechanism that they propose Q , this 
result is consistent with the fact that the magnetic field 
reduces the spatial extent of the electron wavefunctions 
in the directions transverse to the field. This reduces the 
wavefunction overlap of neighboring electrons, resulting 
in a decrease of the electron hopping rate and hence a 
decrease in the rate at which a Coulomb gap forms. This 
is confirmed by our calculations. 

Our model of the Coulomb glass follows that of Bara- 
novskii, Shklovskii, and A.L. Efros (BSE) g. In this 
model, the electrons occupy the sites of a periodic lattice, 
and the number of electrons is half the number of sites. 
Each site has a random onsite energy (j^i chosen from a 
uniform distribution extending from —A to A. Thus, go, 
the density of states without interactions, is flat. A site 
can contain or 1 electron. In order to follow the time 
development of the Coulomb gap, we assume that the 
Coulomb interactions are turned on at time t = 0. The 
Hamiltonian can be written as 

g2 

H = Y1 + H ri,nj9{t) (1) 

where the occupation number Ui equals ^ if site i is oc- 
cupied and —5 if site i is unoccupied, e is the electron 
charge, k is the dielectric constant, and the step function 
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0{t) is for t < and 1 for t > 0. 

The Coulomb gap arises because the stabihty of the 
ground state with respect to single electron hopping from 
an occupied site i to an unoccupied site j requires H 
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At low energies large distances play an important role 
and so we have replaced the sum by an integral over r 
in the exponent. The origin is at site i. n{s) is the oc- 
cupation probability of a site with energy e. T(e', e, r) is 
given by (Q) with r,;j replaced by r, Si replaced by e, and 
Ej replaced by e' . 



where the single-site energy £^=0^-1- J^j 

need to subtract from the density of states those states 

which violate this stability condition. This leads to [p|JlO| 
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F« = l,n' =O)0(i-T,j (£',£, r,,)) 
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e', and is 
= 1 if site i is 



where the single-site energy £j = e 
the lattice constant, n'^ — rii + 1/2; 
occupied and if site i is unoccupied. F{n[,n'j) is the 
probability that donors i and j have occupation numbers 
n[ and n'p respectively, while all other sites have their 
ground state occupation numbers n'^.. t^^ is the number 
of electrons which jump from site i to site j per unit time. 
9{t — Tij) represents the fact that at time t, the primary 
contributions to the change in the density of states will 
be from those hops for which Tij < t 11 1 . In writing eq. 

we assume that these hops together with phonons 
have equilibrated the system as much as is possible at 
time t. The hopping rate r^"^ is given by |^ 

= 7° exp(-^)[l + N{Ai)]F{n^ = 1, 0) (4) 

where a = kub is the effective Bohr radius of a donor, 
and a_B is the usual Bohr radius {as = h^/me'^). We 
will set the mass m equal to the electron mass so that 
as — 0.529A. iV(A^) is the phonon occupation factor 
and reflects the contribution of phonon assisted hopping. 
We are also allowing for spontaneous emission of phonons 
since we are considering a nonequilibrium situation in 
which electrons hop in order to lower their energy. The 
coefficient 7°- is given by [|j 



2e^ 
Sua 



2 r2. 



2hs 



(5) 



where Ei is the deformation potential, s is the speed of 
sound, and d is the mass density. Following BSE, we can 
derive a self-consistent equation for the density of states 
g{e,t): 

g(e,t) = ffoexp |-iy de'g{s\t)J drAwr^ 

F(n(e) = 1, n(e') = Q)e(— +e- e')9(t - rfe', e, r)) } (6) 



Since it is not clear how the stability condition of eq. 
@) can be applied to finite temperatures, we will confine 
our calculations to the case of T = 0. In this case the 
electron occupation factor F{n'^ — l,n'j = 0) — 1 and the 

phonon occupation factor N{Al) = 0. We can solve eq. 
(^) iteratively on the computer. For the first iteration we 
start with g{s', t) = go, and calculate g{e, t). This is then 
used as the input for g{e',t) in the next iteration. Be- 
cause successive iterations converge by alternating above 
and below the correct answer with decreasing amplitude, 
after the first two iterations we use the average of the 
input and output of a given iteration as the input for the 
next iteration. After 11 iterations the typical difference 
between successive iterations is typically less than 1 part 
in 10^. We set the Fermi energy = 0. Because there is 
particle-hole symmetry, we only need to calculate g{e,t) 
for £ < 0. Figure ||a shows the density of states g{e,t) 
as a function of energy at different times, while Figure 
|l|b shows g{e,t) as a function of time at different ener- 
gies. Notice the development of the Coulomb gap occurs 
over many decades in time. The functional form of the 
time dependence of g{£, t) varies with the energy e and 
with go- For example, at the Fermi energy g{fi,t) ^ hit 
for go = 2 X 10^ states/K-A^ and g{fi,t) - t"°"5 for 
go = 6.25 X 10^ states/K-A'^. After an infinite amount 
of time, the density of states at the Fermi energy fj, goes 
to zero and g{e) ^ . For finite times, g(s) ~ |£| in the 
vicinity of the Fermi energy, though there will be thermal 
smearing at finite temperatures. For different values of 
the initial density of states go, we find the same quali- 
tative behavior as a function of time with the depth of 
the dip go — gifi, t) increasing as go increases. Figure ^ 
shows that the width W of the dip increases with go- Ex- 
perimentally the width W increases with the carrier con- 
centration n 1^. This is consistent with our results since 
the noninteracting density of states go increases with n, 
though other parameters such as k may also depend on 
n. The range of widths in Figure ^ is comparable to that 
deduced from experiment Q . 

The temporal development of the Coulomb gap is qual- 
itatively consistent with the experimental observation of 
the long time relaxation of the conductance after the gate 
voltage Vg has been changed. The exact relation between 
the conductance and the density of states is difficult to 
acertain in this case because the system is not in equi- 
librium. However, it is reasonable to assume that the 
conductance reflects the density of states at the Fermi 
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energy. A well known example is Mott's formula for con- 
ductivity a due to variable range hopping H. 



a = Go exp 



(7) 



where To = a/{kBg{n)a^), a is a numerical constant, and 
f/(/i) is the density of states at the Fermi energy. While 
strictly speaking this equilibrium formula does not apply 
to our nonequilibrium situation, we see qualitatively that 
an increase (decrease) in g{fj.) leads to an increase (de- 
crease) in the conductivity. In the experiments, rapidly 
sweeping the gate voltage Vg varies the chemical po- 
tential without allowing time for equilibration. Relating 
the conductance to the density of states means that the 
sweeps over Vg scan the density of states. To get a qual- 
itative feel for this connection, we will use eqn. (^. We 
identify g{^) with g{e) and use To = a/ (kBg{s)a^)- For 
most of the scan the density of states has the linear form 
g[e) = g{So) + a{e — Eo) where £o and the slope a are 
constants. The exponent of 1/4 is appropriate for this 
case. The experiments on indium oxide |^ were done at 
4.11 K, so we set r = 4 K and use the g{e) shown in 
Figure [|. The result is shown in Figure ||. 

Ovadyahu and PoUak noticed that when a magnetic 
field is applied to their indium oxide films, the relax- 
ation rate decreases and the magnetoresistance is posi- 
tive They attributed this to a reduced hopping rate 
resulting from the fact that a polarized spin cannot hop 
onto a site that is already occupied. Our scenario sug- 
gest an additional mechanism since the magnetic field 
reduces the wavefunction overlap in the direction trans- 
verse to the field. The reduced overlap means a lower 
hopping rate and a longer relaxation time. Shklovskii 
and Efros [p"2|-|l4| studied the effect of a magnetic field on 
variable range hopping. Since the average hopping dis- 
tance far exceeds the mean distance R between impuri- 
ties, the hopping electron scatters from many other donor 
sites. As a result, when the magnetic field is transverse 
to the direction of tunneling, the wavefunction decays as 
exp(— p/6). Here we have adopted cylindrical coordinates 
with the magnetic field along the z axis and p is the ra- 
dial coordinate transverse to the z axis. The parameter 
b = A/ 1 In A I where the magnetic length A = y/ch/eH, 
and A describes the scattering and depends on the mag- 
netic field H. For all values of H, Shklovskii has 
shown that b, and hence the wavefunction overlap, de- 
creases as the magnetic field increases. The functional 
form of A depends on the strength of the field. In the 
indium oxide experiments Q = 9 Tcsla. Since this is 
in a weak field regime where A S> a and R <^ A'^/a, we 
can make the approximation b w a[l — (a/A)"*/^] where a 
is the effective Bohr radius. For a weak field we expect 
the hopping rate to go as 7°- exp[~2r^j/ f{z/pij)] 

where the function /(O) = b and f{oo) = a. It is diffi- 
cult to make a quantitative comparison to the density of 



states at = since we do not know the prefactor 7° . 
To get a qualitative feel for the effect of the magnetic 
field, we can use the H = form of 7,°- found in eqn. (||), 
and in the equation (^) for r^"^, we replace exp(— 2rij/a) 
with exp(— 2rij/6). This is a reasonably good approxima- 
tion since for a field of 9 Tesla and a dielectric constant of 
10, a = 5.29 Aand b = 5.162 A. The T = result for the 
density of states at _ff = and H = 9 Tesla is shown in 
the inset of Figure |^. Both curves start at the same value 
oi g{e = ji) — go at time t — Q. Notice that the curve 
at 9 Tesla is slightly above the zero field curve indicating 
that the relaxation is slower in a magnetic field. This 
is qualitatively consistent with experiment jij. To dif- 
ferentiate between our mechanism and spin effects, note 
that our mechanism predicts that the magnetoresistance 
is greater when the field is perpendicular to the current 
than when it is parallel, whereas the magnetoresistance 
should be isotropic in the field if spin effects dominate. 

The experiments found that the dip in the conductance 
as a function of gate voltage Vg and the long relaxation 
times of the conductance following a change in Vg were 
present only at low temperature. These features were 
not observed at T ~ 20 K. From eqs. (||) and (||), we 
see that an increase in temperature will affect g{e,t) in 
two ways. First the thermal smearing of the occupa- 
tion factor F{n[, n'^) will fill in the Coulomb gap to some 
extent. Secondly, as the number of phonons increases 
with temperature, there is an increase in the phonon as- 
sisted hopping of electrons. We expect that this leads 
to a rapid rearrangement of electrons on time scales that 
are too short to observe experimentally. As a result, no 
dips in the conductance and no long time relaxation were 
seen experimentally at higher temperatures. It is difficult 
to calculate these effects because it is not clear how to 
generalize the stability condition (|^) to finite tempera- 
tures. In addition the system is not in equilibrium and 
hence temperature is not well defined for the electrons. 
However the absence of the conductance dips at higher 
temperatures is consistent with our scenario. 

To summarize, we have shown that the time develop- 
ment of the Coulomb gap in a Coulomb glass can in- 
volve very long time scales due to electron hopping and 
rearrangement. These results are consistent with con- 
ductance experiments on disordered semiconducting and 
metallic films. Although we have only considered single 
electron hops, these hops are dependent upon previous 
hops of other electrons through their cumulative effect 
on the single particle density of states. We expect multi- 
electron processes to also contribute to the conductance, 
particularly at long time scales. 
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FIG. 1. (a) Density of states g{e) as a function of energy 
for different times, (b) Density of states as a function of 
time for various energies. Parameters used are Qo = 2 x 10^ 
states/K-A^ T = 0, A = K, k = 10, d = 7.18 g/cm^ 
s = 5.0 X 10^ cm/sec, Ei = 5x 10^ K, and = 4 A. The den- 
sity d is chosen to be that of In2 03. The energy is measured 
from the Fermi energy jj, — 0. 
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FIG. 2. (a) Dimensionless conductivity a joo as a function 
of energy for different times, (b) Dimensionless conductivity 
a joo as a function of time for various energies. The conduc- 
tivity of both (a) and (b) are calculated using Mott's formula 
(|) with T = 4 K, a = 5.29177A and a = 2.23. The rest of 
the parameters are the same as in Figure El 



FIG. 3. Width of the Coulomb gap as a function of go at 
t — 10'''* sec and at infinite time. The width is measured 
at halfway between the minimum and the maximum go 
in the density of states. Inset: Density of states gie) as a 
function of time for different values of the magnetic field H. 
For both curves the energy e = \i and go = 6.25 x 10~^ 
states/K-A^. The rest of the parameters are the same as in 
Figure [ll for both the inset and the main figure. 
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